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Abstract 



> 

[ The neutralino-nucleon (x°-iV) scattering is an important process for direct dark 

matter searches. In this paper we discuss one-loop contributions to the cross section 
■ in the wino-like and Higgsino-like LSP cases. The neutralino-nucleon scattering 

t^J- . mediated by the Higgs x°X° an d Zx°X° couplings at tree level is suppressed by 

the gaugino-Higgsino mixing at tree level when the neutralino is close to a weak 
eigenstate. The one-loop contribution to the cross section, generated by the gauge 
£pH, interaction, is not suppressed by any SUSY particle mass or mixing in the wino- 

Q_i! and Higgsino-like LSP cases. It may significantly alter the total cross section when 

a^o N ~ 10" 45 cm 2 or less. 

> 

X 
S3 



1 Introduction 



Dark matter mass density in the Universe is now measured very precisely by cosmological 
observations, Qm = 0.27±0.04 [1][2]. Now one of the important questions regarding to the 
dark matter is the constituent. The minimal supersymmetric standard model (MSSM) 
predicts the stable lightest supersymmetric particle (LSP) if the R parity is conserved. 
This is one of the attractive features of the model, because the neutralino LSP is a good 
candidate of the dark matter in the Universe because it is a weakly-interacting massive 
particle (WIMP). 

The neutralino LSP is a linear combination of gauginos (bino and wino) and Higgsi- 
nos, which are superpartners of the gauge and Higgs bosons in the SM, respectively. The 
bino-like or Higgsino-like neutralino LSP is the dark matter candidate in the minimal su- 
pergravity model (MSUGRA). Recently many authors have investigated the cosmological 
relic density of the bino-like dark matter in the MSUGRA. It is found that the thermal 
relic density of the LSP is too large compared to the current observations, unless the coan- 
nihilations with other SUSY particles enhance the effective neutralino annihilation cross 
section at the early universe or entropy production reduces the number density after the 
decoupling of the LSP. The anomaly mediated SUSY-breaking model [3, 4] and the string 
models with moduli dominated SUSY breaking [5] predict the wino-like or Higgsino-like 
neutralino LSP. The thermal relic density of the neutralino LSP is too low in the models 
unless the LSP mass (m^o) is above 1 TeV. However, decay of gravitino or other quasi- 
stable particles may produce the dark matter non-thermally so that the relic abundance 
is consistent with the observation [6, 7]. Also, while the LSP with the mass heavier than 
about 1 TeV may lead to the naturalness problem, it may be consistent in the split SUSY 
model [8]. 

Many experiments are now searching for the direct or indirect evidence of the dark mat- 
ter. The counting rates in the direct search experiments depend on the LSP neutralino- 
nucleon (x°-iV) interactions. The cross section above 10~ 42 cm 2 is now explored for 
m^o < 1 TeV. While the annual modulation observed by the DAMA experiment corre- 
sponds to the x°-iV spin-independent cross section around 10~ 42 cm 2 [9], the recent result 
of the CDMSII rejects whole of the DAMA signal region [10] if the spin dependent part of 
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the interaction is negligible. The sensitivity to the dark matter signal may be improved 
up to 10~( 45-46 ) cm 2 or more in future. 

The x°-iV scattering cross section is sensitive to nature of the LSP and the SUSY 
particle mass spectrum. The Higgs and Z boson exchanges are dominant contributions 
to the spin-independent and spin-dependent interactions responsible to the scattering, 
respectively, in the wide parameter region. They are suppressed by the gaugino-Higgsino 
mixing at tree level. When gauginos or Higgsino is much heavier than the weak scale, 
the LSP is close to a pure weak eigenstate and the scattering cross section is strongly 
suppressed. The squark exchange also contributes to the cross section, and it is not 
suppressed by the mixing. However, it tends to be subdominant due to the heavier 
squark masses in the typical models. 

In this paper, we evaluate the one-loop radiative corrections to the wino- and Higgsino- 
like LSP scattering cross sections on a nucleon, which are induced by the gauge interaction. 
These one-loop corrections to the cross section are not suppressed by the Higgsino-gaugino 
mixing. In addition to it, the loop integrals are only suppressed by the weak gauge boson 
masses, because the chargino, which is the SU (2) partner of the LSP, is degenerated with 
the LSP in mass. The gauge-loop correction can dominate the total cross section in a 
limit that the LSP is almost a pure weak eigenstate, setting the "lower limit" of the total 
scattering cross section. When the LSP is bino-like, the one-loop correction is negligible 
since bino does not have gauge charges. 

We note that similar phenomena sometimes appear in radiative correction to the 
SUSY processes. For example, it is known that the mass difference between the LSP and 
chargino may be dominated by the radiative correction due to the gauge loops in the case 
of the wino-like LSP. The tree-level mass difference is 0{m%/Mg USY ). On the other hand, 
the radiative one is not suppressed by any SUSY particle mass, and it is proportional to 
a-ivriw- This is because the loop momentum around m\y dominates in the loop integrals 
due to the mass degeneracy between the LSP and the chargino. In addition, the pair 
annihilation cross sections of the wino- and Higgsino-like neutralino LSPs to two gammas 
at one-loop level are not suppressed by the neutralino mass [11]. It is rather enhanced by 
a non-perturbative effect when the mass is larger than the m-wjai [12]. This effect also 
comes from the mass degeneracy between the LSP and the chargino. 
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This paper is organized as follows. In Section 2, we set up the formula for the general 
low-energy effective action and present the x°-iV scattering cross section. The interactions 
are classified into the spin-dependent and spin-independent ones at the non-relativistic 
limit of the neutralino LSP. We note that there are two classes for the spin-independent 
contributions; one is proportional to the scalar operator of quark, (N\m q qq\N) , and the 
other is proportional to the twist-2 operator, (N\qi(d^ v + <9„7 M — l/2g^ jd)q\N). The 
tree-level contribution to the twist-2 operators, which is induced by the squark exchange, 
is negligible [14]. However, the twist-2 operators is sizable at one-loop level. 

In Section 3, we briefly summarize the tree-level x° couplings responsible to the x°-N 
scattering for the wino- and Higgsino-like LSPs. In Section 4, the gauge-loop correction 
to the effective action in the wino-like LSP case is presented. The dominant contribution 
comes from the W boson loops in this case. We identify the sources of the one-loop 
corrections to the spin-independent interaction, correction to the Higgs coupling of the 
LSP and those to scalar and twist-2 operators. The corrections are only suppressed by 
one-loop factors, not by the Higgsino-gaugino mixing nor any SUSY particle mass. The 
one-loop correction in the Higgsino-like LSP case is summarized in Appendix since the 
structure of the radiative correction is similar to the wino-like case. 

In Section 5, we present some numerical results. Among the gauge-loop contributions, 
sign of the correction to the twist-2 operator is opposite to those to the scalar operator 
and Higgs boson vertex in the case of the wino-like LSP. The cancellation reduces the total 
correction to the spin-independent cross section in wide region of the MSSM parameter 
space. Because of that, the spin-independent cross section induced by the gauge loop 
alone is small, only around ~ io~( 46 ~ 47 ) cm 2 in the limit when the tree-level contribution 
to the cross section is negligible. The total cross section, including the contributions 
at tree and one-loop levels, will be affected by the gauge-loop corrections when the cross 
section is close to the sensitivities of the proposed dark matter search experiments (cr^ojv ~ 
1 0~( 45-46 )cm 2 ). For the Higgsino-like LSP, the cross section induced by the gauge-loop 
diagrams alone is one order of magnitude smaller than that of the wino-like LSP. Section 
6 is devoted to conclusion and discussion. 
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2 Effective Lagrangian for neutralino and nucleon 
scattering 

In this section we present the effective Lagrangian for the x°-iV scattering and the cross 
section formula. The effective interactions of the neutralino LSP with light quarks and 
gluon at the renormalization scale JIq ~ m p are given in a limit of the non-relativistic 
neutralino as follows 

C cS = £ Cf + Cf, (1) 

q=u,d,s 

where 

£f = d Q X°Yl5X° <hV75<? + f q m q x°X° M 

+ ^~ X id^X° 0% + xWWx O^, (2) 

(1) (2) 

+ — ¥id»Yx° o% + %- (<s")x° o%. (3) 

Here, we include terms up to the second derivative of the neutralino field. The first term 
of £^ ff is a spin-dependent interaction, while the other terms in £* s and C c ^ are spin- 
independent "coherent" contributions. The third and fourth terms in Cf and the second 
and third terms in Cf depend on the twist-2 operators (traceless part of the energy 
momentum tensor) for quarks and gluon, 

°lu = \w i^ixlv + dvlp, - ^g^u q, 

0% = {GYG a pu + \g^G a a ,G a ^y (4) 

The scattering cross section of the neutralino with target nuclei is expressed compactly 
by using the coefficients given in and Cf as follows [13], 



A / \ 2 

4 / m^omr \ 
a = - ' 1 



[n p f p + n n f n ) 2 + 4^-^ (a p (S p ) + a n (S n ))' 



J 



(5) 



7i \m^o + rriT / 

The first term in the bracket comes from the spin-independent interactions while the 
second one is generated by the spin-dependent one. 
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In Eq. (5), uit is the target nucleus mass and n p and n n are proton and neutron num- 
bers in the target nucleus, respectively. The spin-independent coupling of the neutralino 
with nucleon, /jv (N — p,n), in Eq. (5) is given as 

f N /m N = £ (/,/T g + ^(?(2,^)+^2,^))(4 1 )+4 2 ))(/io)) 

q=u,d,s 

- ^-JtoSq + \g(2, fit) (gg + gg ] ) (ft), (6) 

where the matrix elements of nucleon are expressed as 

h q = (N\m q qq\N)/m N , 
f TG = 1 — ^2 fr q , 

u,d,s 

lib J\J TI 

(N(p)\0%\N(p)) = J-(p^-I m ^)G(2,/io). (7) 

Here, q(2,p,l), q(2,p,l) and G(2,/2q) are the second moments of the quark, anti-quark and 
gluon distribution functions, respectively, 

?( 2 ,£o) + ?( 2 ^o) = / dx x [q(x,fi 2 ) + q(x,fil)}, 

Jo 

G(2,p*) = I' dx x g(x,j%). (8) 
Jo 

The constant (N — p,n), which is responsible for the spin-dependent contribution, 
is defined as 

a N = d q Aq N , (9) 

q=u,d,s 

2s^Aq N = (N\q W q\N), (10) 

where s M is the nucleon's spin, while (SV) — (A\Sn\A) in Eq. (5) is the expectation value 
of the third component of the spin operator of the proton or neutron group in the nucleus 
A. 

The above formula does not contain heavy quark contributions explicitly. After in- 
tegrating out the Higgs and weak gauge bosons and squarks, the effective interactions 
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of the neutralino with the heavy quarks (Q), which contribute to the spin-independent 
interaction, are 

£q = fQ m Q Px° QQ + /q Qi ftQ 

(1) (2) 

The effective interactions (f f) contribute to £^ ff through the heavy quark loop diagrams. 

It is well-known that the matrix element for ttiqQQ in Eq. (11) can be replaced by 
that for -a s /{\2Tt)G a ^G a ^ due to the trace anomaly. Then, 

(N\m Q QQ\N) = ^f TG m N . (12) 

The second term of Eq. (11) is proportional to Qi 0Q. It would reduce to raqQQx ^ ■, if 
the equation of motion for the heavy quark could be applied. However, this is not justified 
because the heavy quark loop diagram induced by the interaction, which contributes to 
the operator x°x°GG, has an UV divergence. Indeed it is found in Ref. [14], where the 
squark exchange contribution to the elastic scattering is evaluated, that an estimation 
of that operator using the equation of motion for the heavy quarks disagrees with the 
explicit full one-loop calculation by a factor of 2. We need to calculate the vertex for 
X°X°GG, which is induced by the heavy quark loop diagrams, in the original theory. In 
this paper we parameterize the matrix element of Qi 0Q as 

(N\Qi0Q\N) = a eS ^f TG m N . (13) 

by introducing a phenomenological parameter a cS . The precise determination of a e fj 
requires evaluation of the higher-order loop diagrams, and it is out of scope of this paper. 

Eqs. (2), (3) and (11) contain the traceless parts of the energy momentum tensor, 
O 9 , and O 9 , whose matrix elements are scale-dependent. The second moments 
of the quark and anti-quark distribution functions, q(2,fi 2 ) and g(2,/i 2 ), are mixed with 
that of the gluon distribution function, G(2,fi 2 ), once the QCD radiative corrections are 
included. Their scale dependences are compensated by those of the coefficients g^ 1 ' 2 ^ and 
gc' 2 \ so that the total cross section is scale-independent. Thus, 



E (QM) + q(2,p*))(gV> +gW)(fio) + G(2,^)(^ 1) + <?g ) )(A)) 

q=u,d,s 
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For proton 


fru 


0.023 


frd 


0.034 


Its 


0.14 


For neutron 


hu 


0.019 


frd 


0.041 


fr.s 


0.14 



Spin fraction 


Au 


0.77 


Ad 


-0.49 


As 


-0.15 



Second moment at fi = mz 


G{2) 


0.48 






«(2) 


0.22 


«(2) 


0.034 


d(2) 


0.11 


d(2) 


0.036 


5(2) 


0.026 


5(2) 


0.026 


c(2) 


0.019 


5(2) 


0.019 


6(2) 


0.012 


5(2) 


0.012 



Table 1: Parameters for quark and gluon matrix elements used in this paper, (i = 
u,d,s) is taken from the estimation in Refs. [16, 17, 18, 19]. The second moments for 
gluon and quarks at fi = mz and the spin fraction are for proton. Those for neutron 
are given by exchange of up and down quarks in the tables. The second moments are 
calculated using the CTEQ parton distribution [15]. 



m q <fj, 



(14) 



In this paper we use the second moments for gluon and quark distribution functions at 
p, = mz, which are derived by the CTEQ parton distribution [15]. 

In Table 1, we show the parameters for the matrix elements, used in this paper. The 
second moments for the up and down quark distribution functions are sizable. We will see 
in Section 3, the term proportional to g(2, fi 2 ) is generated by the one-loop box correction 
involving the W boson exchanges. The f Tq and f T c represent the fractions of the trace 
part of the energy momentum tensor, or the quark and gluino contributions to the nucleon 
mass, as defined as above. The strange quark contribution to the spin-independent cross 
section is dominant among those for light quarks. Also, the term proportional to Jtg m 
Eq. (6) through heavy quark loops is not negligible. 
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3 LSP elastic scattering induced by tree-level LSP 
couplings 



The mass matrix of the the neutralinos and charginos are given by 



Mm = 



( M x 


—mzsw cos (3 
V rn z s w sin P 





M 2 
mzcw cos P 
—m z c w sin (3 



—m z sw cos (3 
mzcw cos f3 


-n 



m z s w sm.(3 \ 
—mzcw sin (3 







/ 



(15) 



M c = 



M 2 \f2m-w sin {3 

\^2mw cos (3 ii 



(16) 



which is written by the (B, W°, Hf, ff 2 ) bases and (W + , H + ) bases, respectively. 

When the lightest neutralino is wino-like (M 2 <C fx, Mi) or Higgsino-like (/i<Mi, M 2 ), 
the mass eigenvalues of the lightest neutralino and chargino are close to each other. They 
are given as 



m 



w 



for the wino case, and 

= n + 
= V- 



M 2 + 



2 "A* 



m|(l + sin 2/5) 
2(/x-M 1 )(/x-M 2 ) 



-(M 2 + /i sin 2/3) + ... 
(M 2 + n sin 2/3) + ... 



;i7) 



(// - M lC ^ - M 2 s^) + ... 



M 9 2 



;(H + M 2 sin 2/3) + ... 



(18) 



for the Higgsino case (// > 0). The lightest neutralino mass eigenstate becomes very close 
to the pure wino or Higgsino in the limit, therefore m^o ~ M 2 (wino like) or ~ ii (Higgsino 
like). The mass difference between the LSP and the lighter chargino 5 C = (m^- —m^o)/m^o 
is also very small if — M 2 \ ^> rnz- The LSP and the lighter chargino form an SU(2) 
triplet state in the wino-like LSP case or vector-like doublets in the Higgsino-like one with 
the second-lightest neutralino, respectively. 

In this section, we discuss the tree-level contributions to the effective interactions in 
Eqs. (2), (3) and (11) in the wino-like and Higgsino-like LSP cases. The x°-iV spin- 
independent interactions are generated by the t-channel exchange of one Higgs boson and 
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the s-channel exchange of one squark, and the x°-N axial-vector interaction is induced 
by the t-channel exchange of one Z boson, respectively. The interactions of the neutralino 
LSP with the Z and Higgs bosons are suppressed by the mixing of gauginos and Higgsino. 
The squark exchange also generates the twist-2 interaction. The twist-2 coupling is not 
suppressed even in a pure gaugino limit, however, the term is proportional to m^ 4 in the 
amplitude. Overall, if the other SUSY particles are much heavier than the LSP and the 
mixing of gaugino and Higgsino is small, the elastic scattering is suppressed, as mentioned 
in Introduction. The one-loop contributions to the x°-iV scattering will be discussed in 
the next section. 

3.1 Spin-independent interaction 

The spin-independent interaction for the x°-iV scattering arises from one Higgs boson or 
squark exchange at tree level. Since the squark contribution is typically sub-dominant, 
we concentrate on the Higgs boson contribution. 

The neutralino coupling with a quark q, f q , from the Higgs boson exchange is given as 

f [Hi — #2 ( C hxx C hqq , c Hxx C Hqq \ / 1Q \ 

m ~ ^l^r + ^irj' (19) 

where 

sin a cos a: 

Chdd = 77, c Hdd = (20) 

cos p cos p 



for down-type quarks and 



cos a sin a: 

Chuu '■ 77 ? CHuu ~. 7) \ ) 

sin p sin p 



for up-type quarks. Here, a and f3 are the mixing angle of the neutral Higgs bosons 
and the vacuum mixing angle, and m h o and m H o are the light and heavy CP-even Higgs 
boson masses, respectively. The tree-level coupling constants of the neutralino LSP with 
the Higgs bosons, Ch xx an d 

°hxx = [(°n)i2 - (0 N y u t w ] [- sin a(0 N )* 3 - cos a(0 N )\ 4 ] , 

°h X x = [(On) i2 ~ (0 N y n t w ] [cosa(OAr)* 3 - sina(Ojv)* 4 ] , (22) 
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where (On) is the neutralino mixing matrix. Here tw = tan9w, cw = cos 9w and 
sw = sm9w with 6 W the Weinberg angle. 

For the wino-like LSP, the couplings and ch xx are gi ven as 

Here, we assume — M 2 \ ^> mz and cos a ~ sin/5 and sin a ~ —cos/?. The latter 
corresponds to a limit of the heavy pseudoscalar Higgs boson mass (m^). The coupling 
constants are suppressed when 3> M z , as expected. When tan (3 is large, the heavy 
Higgs boson contribution may be dominant. While the LSP coupling with the light Higgs 
boson is suppressed by ~ m x omw/n 2 , the coupling of strange quark to the heavy Higgs 
boson is enhanced proportional to tan (3. 
For the Higgsino-like LSP, we get 

9 1 mw n „ 1 m w _ . . 

C "« = t »'2M^R COS2 ' 3+ 2M^R C ° S2 ' 3 - (24) 
/i > (/i < 0). In this paper we take the LSP mass positive by an axial rotation of the 
LSP field. While the couplings are suppressed by the gaugino masses, the suppression is 
moderate compared with the wino-like LSP. 

Assuming that the dominant contribution comes from the light Higgs boson exchange, 
the cross section for the spin-independent x°~P scattering is approximately given as fol- 
lows, 

for the wino-like LSP, and 

for the Higgsino-like LSP. Here we assume M-i <C // (/i <C M% = Mi) for the wino-like 
(Higgsino-like) neutralino LSP. Note that the cross section for the the wino-like LSP is 
suppressed by /i~ 4 . For /j, — 1 TeV and M 2 = 100 GeV and tan/5 = 10, the cross 
section reduces down to 10~ 46 to 10" 47 cm 2 . It is found in Ref. [20] that the wino- and 
Higgsino-like LSP masses are about (1-2) TeV and the spin-independent cross sections 
are 10~^ 44_48; cm 2 , when imposing the thermal relic density constraint. 
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3.2 Spin-dependent interaction 

The spin-dependent interaction of the LSP arises from one Z boson or one squark ex- 
change. We ignore the squark contribution here, again. The tree-level contribution to d q 
in Eq. (2) from the Z boson exchange is represented as 

d « = si^* (27) 

where T| is for the isospin of a quark q. The tree-level LSP coupling to Z boson, cz xx , * s 

c Zxx = \(0 N ) 13 \ 2 -\(0 N ) U \ 2 . (28) 
For the wino-like LSP, cz xx becomes 

cz xx - cos 2(3. (29) 

On the other hand, in a limit of the Higgsino-like LSP, 

for /j, > (/j, < 0). Thus, again, the coupling for the wino-like LSP is more suppressed 
than that for the Higgsino-like one. 

Using Eqs. (29) and (30), the spin-dependent cross section with proton is approxi- 
mately given as 

o-sd ~ 2 x 10- 38 cm 2 x (Yookev) ' ^ 2(3 (31) 
for the wino-like LSP, and 

o- SD ~ 8 x l(T 39 cm 2 x ( M * ) 2 ( V- ) " 2 cos 2 2(3 (32) 

VIOOGeV/ VIOOGeV/ 1 V ; 

for the Higgsino-like LSP. Here, we take Mi = M 2 . The analysis in Ref. [20] shows that 
the spin-dependent cross sections are 10~^ 41_45 ^cm 2 for the wino- and Higgsino-like LSPs 
when imposing the thermal relic density constraint. 
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Figure 1: (a) Diagram contributing to the anomalous Higgs boson vertices of the neu- 
tralino and (b) box diagram contributing to the x°-iV scattering in the case of the wino- 
like LSP. 

4 Elastic scattering induced by one-loop effective ac- 
tion 

In the previous section, we discuss that the interactions responsible for the x°-N scattering 
are suppressed by the gaugino-Higgsino mixing or squark masses at tree level. However, 
this is not true for the radiative corrections to the effective interactions if the dark matter 
is wino- or Higgsino-like, because of the mass degeneracy between the LSP and its 577(2) 
partner. In this section, we derive radiative corrections to the effective interactions in 
Eqs. (2), (3) and (11), and it is found that some of them are only suppressed by the weak 
gauge boson mass at most. 

We first discuss the anomalous Higgs boson vertices of the neutralino and the box 
diagram contributions involving the W bosons to the effective interactions for the case of 
the wino-like LSP. The numerical result will be shown in the next section. For the case 
of the Higgsino-like LSP, we present the explicit formula for the radiative corrections in 
Appendix. 

The gauge interactions of the wino-like neutralino and chargino are 

Ant = -— (F7 M r^+^.)+e— Fyr^+eryr^. (33) 

Sw y ' Sw 

Here we ignore the mixings of the neutralinos and charginos for simplicity. These interac- 



12 



tions induce the anomalous Higgs boson vertices and the box diagrams which contribute 
to the x°-iV scattering. The Feynman diagrams are given in Fig. 1. The radiative cor- 
rection to the Z boson vertex is not induced at the one-loop level due to the Majorana 
nature of the LSP. 

We start from the radiative correction to the the Higgs boson vertices. The tree-level 
contribution is given in Eq. (22). The radiative corrections to Ch xx an d c h X x a ^ one-loop 
level are expressed as 

fchix = ^Ma-^lF^M + ScF^M], (34) 
Sc H X x = -^cos(a-/3) F^\x w ) + 6 c f£\x w ) , (35) 
respectively. The correction to f q [H](— f* ree [H] + 8f q [H\) is therefore 

5f q [H] = -f- ( + 5cH *f H A . (36) 



Here, xw = rn^jm 1 ^ and b~w = \f^~ x w /A. We expand the radiative corrections by 
5c = (jn x - ~ m^o)/m^o. When the gaugino and Higgsino masses are comparable to 
the weak gauge boson masses, this expansion of 5c is invalid, however, the one-loop 
contributions to the cross section are negligible. We are interested in a case where the 
gaugino or Higgsino mass is much larger than the weak gauge boson masses and the LSP 
is close to the Higgsino or wino weak eigenstate. This expansion is justified in the case. 
The mass functions are 

F$\x) = ^-(2 + x(2-x))t a n- 1 &)-2^(2-x\og(x)), (37) 
b w ^x 



F<« W = |- tl m-'(^|)-^-L(2 + ,). (38) 

Since F^\x) becomes 2it in a limit of x — > 0, the correction does not vanish in the heavy 
neutralino limit. When the pseudoscalar Higgs boson mass is large, 5ch xx is maximum 
and 5chxx i s vanishing. 

The effective interactions in Eqs. (2) and (11) receive the radiative correction from 
the the box diagrams (Fig. 1(b)) at the renormalization scale fi ~ mw of the form, 

5£f[box] = H[box] F7 M 75X° 97,759 

+ 5/Jbox] m q x°X° m + <5/o[box] x°X° ¥ ^ 
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^Yx° 01 + ^^F(^)(^)x° 01, (39) 

where 



Oi Ot Ot 

<WJbox] = — f-F AV (a; w ), 5/Jbox] = — |_F S i(x w ), 5/^ [box] = -^-F S2 (x w ), (40) 

772^/ m W 

^(i)[ b0 x] = Sg^[box] = *L Ft2 (x w ). (41) 

In our calculation of the radiative corrections, we ignore 0{m 2 ), and expand the loop 
integrals up to order of p and 5c, where p is the quark external momentum. In this ap- 
proximation the loop integrals are expressed analytically by B functions and its derivatives 
[21]. This procedure is not justified for box diagrams with the external or internal top 
quarks since it is heavier than the weak gauge bosons, however, the radiative corrections 
should be suppressed by the top quark mass. Thus, they are sub-dominant compared 
with the other lighter quark ones. 

The loop function can be expanded as Fj(x) = Fj°\x) + 5cFj 1 \x) (I = AV, SI, S2, 
Tl, and T2), and they are given as follows; 

F§(x) = ^^(8-x-x 2 )t a n-\^)-^x(2-(3 + x)\og(x)), 

F av( x ) = Tl^rv^tan-^ 21 '' ' 



Ab 3 w v v v^ 7 2b 2 w ' 

F£\x) = 4?(x) = 0, 

^s ( 2 0) (*) = -^(2 + a?)tuT\^)-±,/Z{l-2x-x(2-x)]Dg(z)), 
24 \ x 9o 



+-i^(2(6 -x)+6 log(25 c ) - (3 - x 2 ) log(x)), 

24 yJX 

F$(x) = h w {2 + x 2 )t^ 1 &) + ^- y ^(l-2x-x(2-x)\ogx), 

F$(x) = -^(l-2x + x*)t^&) 

+ 7 ^ r (2x + 6 log(25 c ) - (3 + x 2 ) log(z)), 
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Figure 2: A two- loop diagram contributing to GGx°X°- 



4b w 
1 



246 3 w 



x(2 — Ax + x 2 ) tan 1 
(16 + 30x- 30x 2 + 



,26 



ir 



— — 2x — x(2 — x) log(a;)), 



5x 3 ) tan 1 



.26 



X 



246^ 



a/x(28 - 10a; - 5x(4 - x) log(x)). 



(42) 



Note that F$\(x) is zero due to the chiral nature of Wqq vertex. For the Higgsino-like 
LSP, it does not vanish because the Z boson couples with both qi and qn. See Appendix. 

The functions F^ (x) and F^i (x) are non- vanishing even if x approaches to (or m^o 
is increased), and they become — 7r/24 and 7r/6, respectively. Thus, in addition to the 
correction to the Higgs boson vertices in Eq. (36), those to the scalar and twist 2 operators 
induced by the W boson loops do not vanish in a heavy LSP limit, and they contribute to 
the spin-independent cross section. On the other hand, the spin-dependent cross section 
depends on Fjw(x), which is suppressed in the heavy LSP limit as cx mw/m^o. 



The functions, F, 



x) 



Fg\x) 



and F$(x) 



are proportional to 1/y/x for small x. 



This does not cause a problem. The mass difference 5c is proportional to x% 2 at tree 

1/2 

level, and it becomes proportional to a-i^w due to the radiative correction when the LSP 
mass is much larger than the W boson mass. Thus, Sc/y/xw ~ %w or ~ a 2 - Therefore 
the perturbation by 5c is not broken. 
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5 Numerical results 



In the previous sections, we derived the formulas for the one-loop corrections to the 
effective interactions. We now calculate the cross section following the procedure in 
Section 2. 

First, we discuss the cross section for the spin-independent x°~P scattering induced 
by the gauge-loop diagrams presented in the previous section, assuming the tree-level 
contribution to the cross section is negligible. The total cross section involving both the 
tree and one-loop level contributions is more sensitive to the MSSM model parameters, 
and it is discussed later. 

For the wino-like LSP, the spin-independent interaction induced by the W boson at 
one-loop level is approximately represented by 

x , na%m N ( f T f TG \ na 2 2 m N ( f T f TG \ 

SfN = IJL ~ + ~ j " "4T 1-5, 24 + ^324 J 

+ E !i ^(9(2,n4)+9(2, n 4)). (43) 

q=u,d,s,c W 

Here, we take the pseudoscalar Higgs boson and the LSP to be much heavier than the W 
boson mass. The first term comes from the correction to the light Higgs boson vertex, 
proportional to 5ch xx - The second one is from those to the scalar operators and it is given 
by Sf q and Sf' q in Eq. (39). The third (last) one represents those to the twist-2 operators, 
which are proportional to 8g q ^ and Sgj® in Eq. (39). Here, we include the top and bottom 
quark contributions to the Higgs vertex correction, while those to the scalar and twist-2 
operators, which are suppressed by the top quark mass, are ignored. As discussed before, 
a cS is a phenomenological parameter for (N\Qi flQ\N). Eq. (43) is numerically given as 
(in the order of the terms in Eq. (43)) 

5f N = -5.8 x 1(T 10 - (5.3 + 1.6aeff) x KT 11 + 3.9 x 10" 10 (GeV" 2 ) 

= -2.6 x l(T 10 (GeV~ 2 ), (a cS = 1). (44) 

Here, m h o = 115 GeV and the parameters in Table 1 for the hadronic matrix elements 
are used. We assume a cf j = 1 in the second line. We found the contributions from the 
terms with different spin structure cancel each other. This value for 5 fx corresponds to 
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<7si — 3.0 x 10~ 47 cm 2 , if the tree-level contribution is negligible. The correction is still 
large enough to alter the total cross section close to the proposed sensitivities for the 
future experiments (lcH 45-46 ) cm 2 ). Here, we present the spin-independent interaction of 
the neutralino with proton, however, the value for that with neutron is almost the same 
as it since the numerical difference comes only from the Higgs exchange contributions. 

In Fig. 3, we show the one-loop induced cross section for the spin- independent x°-p 
scattering as a function of fas for a c g = 1.0 and 2.0. The tree-level contribution is assumed 
to be negligible. The parameters fr s and a c g have large theoretical uncertainties. The 
precise determination of a e fr requires calculation of two- loop diagrams, such as in Fig. 2. 
Here, we take m^o = 1600GeV and the parameters in Table. 1 except for f Ts . Also, we 
use FeynHiggs [22] in order to calculate the Higgs boson masses and mixing angle. In 
this figure, we use tan/5 = 10, m j 4=1000GeV, m stop =2000GeV, and A top —0, which lead 
to m h o = 116GeV. The cross section is very sensitive to f Ts since the correction to the 
Higgs boson vertex is one of the dominant corrections. When f Ts is smaller than 0.1, the 
cross section is significantly suppressed. Since the box diagram corrections to the scalar 
operators are numerically small, the a e g dependence is relatively small. 

In Fig. 4 the one-loop induced spin-independent cross section is presented as a function 
of tan/5 and vtia- Here, we take m^o=1600GeV, m stop =2000GeV, A top =0, a c s = 1, and 
the parameter in Table. 1. The tan/5 dependence is moderate. The one- loop induced 
cross section rises as large as 10~ 46 cm 2 when the pseudoscalar Higgs boson is light. 

In the case of the Higgsino-like LSP, the one-loop induced cross section for the spin- 
independent scattering is smaller than that in the wino-like case since the SU(2) charge 
is smaller. When the pseudoscalar Higgs boson and the LSP are much heavier than the 
W boson mass, S/n is given as 

5f N = 4.3 x 10~ n - (5.1 + 0.1a cff ) x 10" 11 + 6.0 x 10~ n 

= 5.1 xlO- 11 , (a cff = l). (45) 

Here, we take m^o = 115GeV, again. The order of the terms in the first line is the same 
as Eq. (44). This value for Sf N corresponds to <r SI ~ 1.1 x 10~ 48 cm 2 when the tree-level 
contribution is negligible. 

Now we discuss the gauge-loop correction in the total cross section in the general 
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12 



CM 




0.05 0.1 0.15 0.2 0.25 0.3 

fls 

Figure 3: One- loop induced cross section for spin-independent x°~P scattering as a func- 
tion of frs for a c g = 1.0 and 2.0 in the case of the wino-like LSP. We assume that the 
tree-level contribution is negligible and the process is induced by the one-loop diagrams. 
Here, we take m^o = 1600GeV and the parameter in Table. 1 except for Also, we 
use tan (3 = 10, m/i=1000GeV, m stO p=2000GeV, and A top =0, which determine the light 
Higgs mass and the mixing angle. 
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Figure 4: One-loop induced cross section for spin-independent x°-p scattering is pre- 
sented as a function of tan/3 and uia- The wino-like LSP is assumed. Here, we take 
m^o = 1600GeV, m stop =2000GeV, A top —0, a c s = 1, and the parameter in Table. 1. 

MSSM parameter space. In above paragraphs we showed the cross section induced by the 
one-loop diagrams alone, assuming the tree-level cross section is negligible. When the tree- 
level amplitude dominates, the one-loop correction to the cross section is approximately 
expressed as ~ 2y / o"i_i oop /cr tr ee, where o" tree and ai-i oop are the tree-level and one-loop 
induced cross sections, respectively. When a"i-ioop/°"trcc ~ 1/100, the correction to the 
total cross section is about 20%. 

In Fig. 5 we show the cross section involving the tree-level and one-loop contributions, 
o"totah and (T to tai/ctree f° r the spin-independent x°~P scattering as as functions of \i and 
Ma- Here, we assume the wino-like LSP and take m^o ~ 200GeV. The upper two figures 
are for tan f3 — 4 and the lower ones are for tan f3 — 40. See the caption for the other input 
parameters. It is found that the radiative correction is about 50% in the plots, where o" total 
is above 10~ 45 cm 2 . The radiative correction is relatively significant in larger tan/3, since 
the coupling of the LSP with the light Higgs boson at tree level is more suppressed in 
the case as discussed in Section 3.1. In Fig. 6, we show the case where m^o ~ 2 TeV. 
It is found that <7 to tai is less than 10~ 45 cm 2 , however it is dominated by the gauge- loop 
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Figure 5: Total cross section for spin-independent x°~P scattering, fJ to tab (left) for the 
wino-like LSP in the (/i, tua) plane. The cross section is depicted in the unit of lCT 47 cm 2 . 
Here, M 2 = 200 GeV, M 1 /M 2 = llt 2 w , m- t = 2 TeV and tan/3 = 4(40) for upper (lower) 
plot. rJtotai includes the gauge-loop contributions while the squark exchange contribution 
at tree level is ignored. In the right, we give the ratio o" to tai/ 'tree for the same parameters. 
(j trcc is the cross section at tree level. 



correction in the plots, since the LSP is very close to the pure wino state. 

In Figs. 7 and 8, we show the effect of the radiative correction for the Higgsino-like 
LSP. The radiative correction is less than ~ 10% for m^o ~ 200 GeV even when Ma 
and Mi are as heavy as 2 TeV. Compared with the wino-like case, the suppression of the 
tree-level Higgs boson vertex by the gaugino-Higgsino mixing is more moderate while the 
radiative correction by the gauge-loop diagrams is smaller. The correction is negative 
(positive) for /i > 0(< 0), relatively to the tree-level contribution. For m^o ~ 2 TeV, the 
one-loop contribution cancels the tree-level contribution to reduce the total cross section 
less than 10~ 50 cm 2 in the figure. 

Finally, we discuss the spin-dependent cross section. As discussed in the previous sec- 
tion, the one-loop contribution to the process is suppressed by mw/m^o in the amplitude, 
contrary to the spin-independent cross section. In Fig. 9 the one-loop induced cross sec- 
tion for the spin-dependent x°-p scattering is given as a function of the LSP mass under 
the assumption of the wino-like LSP. In this figure we assume again that the tree-level 
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Figure 6: Same as Fig.5, but M 2 = 2 TeV, m~ t = 20 TeV. 
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Figure 7: Total cross section for spin-independent x°~P scattering, o- tota i (left) for the 
Higgsino-like LSP in the (/z, m^) plane. The cross section is depicted in the unit of 
l(r 47 cm 2 . Here, = 200 GeV. M 1 /M 2 = 5t 2 w /3, m~ t = 2 TeV and tan/5 = 4(40) for 
upper (lower) figure respectively. In the right, we give the ratio o" tota i/(T tr ee. 
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Figure 8: Same as Fig. 7, but (i = 2 TeV and m~ t = 20 TeV. 



contribution is negligible. The asymptotic behavior for the cross section is 

ffs ° - ^ US/'] 

K UxlO-rfx^)-'. (46) 

Thus, from Eq. (31), it is found that the one-loop correction becomes significant for 
m^o ~ 100 GeV when 1 TeV. For 3> M 2 , the one-loop contribution is constructive 
to the tree-level one for the spin-dependent x°-p scattering while it is deconstructive for 
the scattering. 

For the Higgsino-like LSP, the one-loop induced spin-dependent cross section is one 
order of magnitude smaller than that for the wino-like LSP. It behaves as 

<*■> - L5xlr4Wx (lsSv)" 2 < 47 > 

in a large m^o limit. The tree-level contribution does not suffer from significant sup- 
pression by the Higgsino-gaugino mixing as in Eq. (32), and the one-loop correction is 
negligible as far as the gaugino masses are smaller than 10 TeV. For M 1; M 2 3> 
the one-loop contribution is deconstructive (constructive) to the tree-level one for the 
spin-dependent x°~P (x°~ n ) scattering. 
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Figure 9: One- loop induced cross section for the spin-dependent x°~P scattering in the 
case of the wino-like neutralino DM. Here, we assume that the tree-level contribution is 
negligible. 

6 Conclusion and discussion 

In this paper, we studied one-loop correction to the neutralino-nucleon scattering processes 
for the wino- and Higgsino-like LSPs. The x°-N scattering is relevant to the detection rate 
of the neutralino dark matter. The scattering occurs dominantly through the exchange of 
the Higgs or Z boson at tree level, which is suppressed by the gaugino-Higgsino mixing. 
This is especially the case when the gauginos or Higgsino mass is large. 

On the other hand, the scattering cross section receives the one-loop contribution from 
the processes involving the weak boson exchange. These processes are not suppressed 
neither by the small mixing angle nor by heavy SUSY mass scale, therefore it could be 
significant part of the one-loop corrected cross section in the limits of the wino- and 
Higgsino-like LSPs. The spin-independent cross section for the wino-like (Higgsino-like) 
LSP receives the sizable one-loop correction, when the Higgsino (gaugino) mass is heavier 
than about ITeV and the spin-independent cross section is smaller than about 1CT 45 cm 2 
(10- 46 cm 2 ). 

There has been impressive progress on experimental techniques for the dark matter 
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searches. The on-going experiments are now sensitive enough to exclude some of the 
MSSM parameters which may not be accessible otherwise. Many proposals aim to push 
the sensitivities further. The better determination on the x°-iV cross section is also needed 
to understand the local dark matter density and its velocity distribution. Although we 
calculated only a part of the correction of the scattering cross section, further investigation 
on the loop effect might be needed in future. 

Acknowledgment 

This work is supported in part by the Grant-in-Aid for Science Research, Ministry of 
Education, Science and Culture, Japan (No. 15540255, No. 13135207 and No. 14046225 for 
JH and No. 14540260 and No. 14046210 for MMN). 



Appendix: Higgsino-like LSP 



In this Appendix, we present the radiative corrections to the effective interactions for 
the LSP scattering with nucleon when the LSP is Higgsino-like. The Higgsino-like LSP 
accompanies the chargino and the second-lightest neutralino, whose masses are degenerate 
with that of the LSP. Furthermore, the Higgsino-like LSP has an interaction with the Z 
boson. 

The one-loop corrections to the LSP couplings with the Higgs bosons are given as 



1 «2 



^sin(a-/3) [F^(x w ) + S c F^(x w ) 



1 a 2 



4 4vrc^ 
1 a 2 



sm(a-P) Fg\x z ) + 6 N Fg>(x z ) 



.ir£cos(a-/3) F§\x w ) + S c F^(x w ) 

4 47T L 

+ -t^V cos(a - (3) \f£\x z ) + 5 N Fg\xz) 



*(i), 



4 47TC 2 



w 
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(49) 



The contributions from the box diagrams to the coefficients in the effective Lagrangian 
(39) are approximated as 
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F$\xz) + 5 N F#(x z 



(1 = 1,2). (53) 



-QqS^r for a quark g. The Z boson contributions 



Here, L q = Tg — Q 9 s^ and i? g = 
appear in the above equations. In the equations, xz = m 2 z /m^ , bz = \Jl — ^z/4, and 
5n = ( m x ( 2 ~ m x°)/ m x w ith m^o the heavier neutralino mass. All mass functions are the 
same as in Eqs. (42) except for F^(x) and F^{\x), 



-i-(4 - i(2 - x) j taiT'fe + i^P - zlog(i)), 
' (2-x)tan- 1 (^I) 



(54) 



4fo| 



\by \/.v 



(4 - 4x - 2(4 - x) log(25jv) + (4 - x) log(x)). 



(55) 
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